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^7^ ■ We present the current results of our simulations of AA = 1 supersymmetric 

^ i Yang-Mills theory on a lattice. The masses of the gluino-glue particle, the a-rj' , 

K> I the a~/o meson, and the scalar glueball are obtained at finer lattice spacing 

H ■ than before, and extrapolations towards vanishing gluino mass are made. The 

- - - calculations employ different levels of stout smearing. The statistical accuracy 

as well as the control of finite size effects and lattice artefacts are better than in 
previous investigations. Taking the statistical and systematic uncertainties into 
account, the extrapolations towards vanishing gluino mass of the masses of the 
fermionic and bosonic states in our present calculations are consistent with the 
formation of degenerate supermultiplets. 



1 Introduction 

J\f = 1 supersyinmetric Yang-Mills theory (SYM) is a theory of gluons, described by non- 
Abelian gauge theory, and their superpartners, the gluinos, which are spin 1/2 Majorana 
fermions in the adjoint representation of the gauge group. The on-shell Lagrangian of SYM 
in Minkowski space is 



£ = tr 



1 i — 7T) — 



(1) 



where Ff^y is the non-Abelian field strength formed out of the gluon fields A^(x), \{x) is the 
gluino field, and D^ denotes the gauge covariant derivative in the adjoint representation. 
The supersymmetry of the theory is broken softly by the gluino mass term. 

In several aspects the model is similar to QCD [T]. The importance of related theories 
in extensions of the Standard Model and the connection to QCD [2] is a main motivation 
for the recent interest in these models. SYM is asymptotically free and is assumed to 
show confinement. Thus gluons and gluinos are not particle states in the physical Hilbert 
space, which instead contains bound states of gluons and gluinos. These bound states must 
form degenerate supermultiplets, if supersymmetry is unbroken. The determination of the 
masses of bound states is a non-perturbative problem. Previous work of our group has been 
dedicated to the calculation of the low-lying masses of bound states in SYM by means of 
numerical simulations on a space-time lattice [SJ |3] . The present article is a continuation of 
our preparatory work, which should be referred to for more details and references [S]. 

On the lattice, supersymmetry is generically broken [B] . In SYM a fine-tuning of the bare 
gluino mass parameter in the continuum limit is enough to approach the symmetries of the 
continuum theory 13 |B]. These symmetries include (spontaneously broken) chiral symmetry 
and supersymmetry. The theoretical prediction of the existence of a supersymmetric chiral 
continuum limit needs to be confronted with the numerical lattice simulations. 

A necessary condition for the restoration of supersymmetry in the continuum limit is the 
degeneracy of fermionic and bosonic masses. From low energy effective theories, predictions 
have been made for two low- lying supermultiplets [9l[T0]. Each multiplet consists of a scalar, 
a pseudoscalar, and a fermionic spin 1/2 particle. The lighter multiplet consists of a O"*""*" 
and a '" glueball and a gluino-glue state. The gluino-glue is an exotic spin 1/2 Majorana 
fermion, which can be created by the operator 

Ogg = Y.a^,iT[F^''\], (2) 

with a^y = ^ [7/x,7iy]- The heavier multiplet is built from the scalar meson a-/o, represented 
by AA, the pseudo-scalar meson a-r/', represented by A75A, and a gluino-glue state. 

In our previous work [3l U] the expected degeneracy of the fermionic and bosonic masses 
was not observed. The gluino-glue appeared heavier than its lightest possible superpartners. 
In that work, however, the masses were obtained at a fixed lattice spacing, and the influence 
of finite volume effects was not known with sufficient accuracy. In order to obtain results 
relevant for the presumed supersymmetric continuum limit, the chiral limit, the continuum 
limit, and the infinite volume limit have to be extrapolated from the simulations. 

In our preparatory work [5] we have made a detailed analysis of finite volume effects, and 
observed supersymmetry breaking effects due to the finite spatial extent of the lattice. The 



results of these calculations allow to estimate the lattice sizes necessary for neglecting them. 
Moreover, simulations at a lattice spacing smaller than before indicate that the influence 
of the finite lattice spacing is larger than expected, and could be the cause of the apparent 
supersymmetry breaking. 

In the present article we conclude our investigations of the finer lattice spacing, taking 
into account also different levels of stout smearing. The masses of the gluino-glue particle, 
the a-r]', the a-/o meson, and the scalar glueball are obtained with the high statistics that 
turned out to be necessary, and extrapolations towards vanishing gluino mass are made. 
The statistical accuracy as well as the control of finite size effects and lattice artefacts is 
better than in all the previous investigations. 

2 Numerical simulations 

In our numerical simulations of SYM, the Euclidean version of the model is formulated on 
a space-time lattice with an action proposed by [7]. We are considering the case of gauge 
group SU(2). The gauge field dynamics is defined by the tree-level Symanzik improved 
plaquette action. The inverse bare gauge coupling in the current simulations was /3 = 1.75. 
The results will be compared to the previous ones at /3 = 1.6 [2111]. The gluinos are described 
by Wilson fermions in the adjoint representation. The Wilson-Dirac operator 

X'^w )x,a,a;y,b,/3 

4 

= Sxy6abSa(3 - ^ X! [(^ ~ lfi)al3{V^{x))abSx+fj.,y + (1 + lfi)al3{V^{x - fJ-))abSx-fi,y (3) 
11=1 

contains stout smeared [11] gauge links V^{x) in the adjoint representation. The hopping 
parameter k is related to the bare gluino mass via k = l/{2mg + 8). The recovery of 
both supersymmetry and chiral \]{1)r symmetry in the continuum limit requires to tune 
the hopping parameter to the point Kc{(3), where the renormalised gluino mass vanishes 
[3 IB]- In practice, this is achieved by monitoring the mass of the unphysical adjoint pion 
(a-7r), which can be defined in a partially quenched setup. The adjoint pion mass, which is 
expected to vanish at Kc, can can be numerically obtained relatively easily. Its correlator is 
defined as the connected contribution of the correlator for the a-r]' meson. The consistency 
of the tuning according to ma-7r with the supersymmetric Ward identities was shown in [3]. 
Further details can be found in . 

The configurations have been obtained by updating with a two-step polynomial hybrid 
Monte Carlo (PHMC) algorithm |12| [3]. Near Kc the occurrence of low eigenvalues of 
the Hermitian Wilson-Dirac operator makes it necessary to introduce correction factors to 
the polynomial approximation in the PHMC algorithm. These have been obtained, when 
necessary, from the correct fermionic contribution of the lowest eigenvalues. In contrast to 
the theory in the continuum, the lattice theory has a (mild) sign problem. The Pfaffian 
obtained by the integration of the Majorana fermions can sometimes have a negative sign 
|13] . When necessary, we included this sign in the reweighting. To reduce the statistical 
errors we have chosen the parameters of our present simulations such that the reweighting 
with correction factors and Pfaffian signs is not relevant for the results. 



On the basis of our investigation of finite volume effects [5] , the lattice sizes have been 
chosen such that finite volume effects can be neglected in comparison with statistical errors. 

For a comparison of dimensionful quantities we use the Sommer parameter ro, obtained 
from the static quark potential. For illustration, QCD units are being used by setting the 
Sommer parameter to ro = 0.5 fm. We have always taken the value of tq obtained by an 
extrapolation to Kc of the data obtained at «; > Kc- In QCD units the lattice spacing is 
a = 0.058 fm at /? = 1.75, and a = 0.088 fm at /3 = 1.6 for the runs with one level of stout 
smearing. A convenient substitute for the gluino mass is the squared adjoint pion mass in 
physical units, (rorria-Tr)^. 

The parameter sets, underlying our calculations at /? = 1.75, are summarised in Tab. [TJ 
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Table 1: Summary of the simulation parameters. Ig is the level of stout smearing in the 
Wilson-Dirac operator. The value of rg/a used in (rorria-Tr)^ is obtained from 
the extrapolation to the chiral limit. Its value is ro/a = 9.02(18) for Ig = 1, 
and ro/a = 8.663(81) for l^ = 3. A'^conf is the number of configurations in the 
measurements of the pion mass. The number of reweighting factors smaller than 
0.98 is around 3% at k = 0.1495 and around 1% at k = 0.1494 (/^ = 1; 323 ^ 54^ _ 
In all other runs the reweighting is not taken into account because the number of 
negative Pfaffians is always below 1%. Tp is the integrated autocorrelation time. 



3 Bound state masses 



We have investigated the masses of the gluino-glue particle, the a-rj' and the a-/o mesons, 
and the scalar glueball. The masses of the bound states are obtained from fits to the 
corresponding correlation functions. In case of the a-vr and the gluino-glue the correlation 
functions yield rather good fits for a number of fit intervals. By means of a histogram 
method [H] reliable mass estimates for these particles could be obtained. In case of the 
mesons a-rj' and a-/o the mass estimate is based on a single optimally selected interval in 



Euclidean time. This interval is fixed to the region where a plateau of the effective mass is 
observed. A second estimate is taken from a fit interval of the same length, but shifted by 
one positive unit in Euclidean time. The difference of these two masses provides an estimate 
of the systematic error of the procedure. Further details of the histogram method and the 
fit procedure are explained in [S]. The fit procedure for the glueball is detailed below. The 
results for the masses are collected in Table |2j 
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Table 2: The masses of the relevant bound states in units of the Sommer scale rg. For 
the mesons a-r]' and a-/o and for the glueball 0^^ the indicated error is given as 
(statistical error) (systematic error of the plateau estimation). The a-vr and gluino- 
glue {gg) mass is obtained with the histogram method and already includes an 
estimate of the systematic error. 



The extrapolation to the chiral limit is obtained from a linear fit of the masses as a 
function of the squared adjoint pion mass. Up to lattice artifacts and finite size effects this 
limit coincides with the supersymmetric limit, where the particle masses should be grouped 
in the predicted multiplets. 

For the lattices with spatial extent L = 24 our statistics is much higher than for L = 32. 
An extent of L = 24 is also sufficiently large in view of finite size effects. Therefore we have 
decided to include only the L = 24 data in the extrapolations. Adding the L = 32 data 
would lead to small changes only. 

The best accuracy can be obtained for the mass of the gluino-glue particle. Its correlator 
has been obtained using a combination of APE and Jacobi smearing. The gluino-glue 
mass is shown in Fig. [T]as a function of the squared adjoint pion mass, together with the 
extrapolation to the chiral limit. In all the figures the error bars include both the statistical 
and systematic errors. 

The correlators for the a-r/' and the a-/o bosons contain connected and disconnected 
contributions. The disconnected parts have been calculated using the stochastic estimator 
method [T2j, taking into account the exact contribution of the 100 lowest eigenmodes of the 
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Figure 1: The gluino-glue mass as a function of the squared mass of the adjoint pion in units 
of the Sommer scale, and the corresponding linear fit. 

even-odd preconditioned Hermitian Wilson-Dirac operator. The disconnected contributions 
are especially significant at smaller adjoint pion masses. The statistical fluctuations in this 
part are larger than in the connected contribution and lead to a bad signal-to-noise ratio 
in the correlators. In Fig. [5] the masses of the mesons are displayed together with the 
extrapolated values at the chiral limit. For comparison, the figures additionally include the 
linear fit of the gluino-glue mass. 

The operator corresponding to the 0"''^ glueball is a combination of gauge links. The 
correlation functions are as usual afflicted by large statistical fluctuations. We are using 
variational smearing methods to improve the signal [IB]. Smearing the underlying gauge 
links leads to a basis of different operators with the same quantum numbers. A combination 
of these operators providing the best overlap with the particle state is estimated numerically. 
From the correlations of the basic operators the correlation matrix C{t) is obtained. The 
optimisation corresponds to the solution of the generalised eigenvalue problem 



C{t)v = X{t,to)C{to)v 



(4) 



for a given fixed to, in our case to = 0. The resulting eigenvalues A(t, 0) as a function of t 
are taken as input for a fit procedure analogous to the case of the mesons a-r]' and a-/o. 
The quality of the signal is improved considerably by the variational method. However, the 
signal is yet not as good as for the gluino-glue, and the systematic error estimate of the 
plateau estimation is taken into account. The results for the glueball masses are shown in 
Fig.S 
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Figure 2: The a-r]' mass and the a-/o mass as functions of the squared mass of the adjoint 
pion in units of the Sommer scale, and the corresponding linear fit. Also shown 
is the fit for the gluino-glue. 
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Figure 3: The glueball mass as a function of the squared mass of the adjoint pion in units 
of the Sommer scale, and the corresponding linear fit. Also show^n is the fit for 
the gluino-glue. 



Whereas in previous work a significant gap between the masses of the gluino-glue and its 
supposed superpartners showed up, the figures presented here display considerably smaller 
differences. To make the comparison quantitative, we confront the present results with those 
of our earlier work in [3], done at a larger lattice spacing. The masses are summarised in 
Table El They have been converted into units of MeV using the QCD scale setting ro = 0.5 
fm. 
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Table 3: Comparison of the bound state masses in units of MeV, extrapolated to vanishing 
gluino mass, at the two values of /3. The lattice spacing at /3 = 1.6 is between 0.088 
and 0.097 fm depending on the level of stout smearing. In the current simulations 
it is approximately 0.055 fm for Ig = 1, and approximately 0.058 fm for Ig = 3. All 
values are obtained using the QCD units by setting tq = 0.5 fm. 



One can observe significant differences between the results at the larger lattice spacing 
(/? = 1.6) and the current lattice spacing (/? = 1.75). The gluino-glue mass gets smaller 
and the other masses larger when the lattice spacing is reduced. Taking the statistical and 
systematic uncertainties into account, the extrapolations towards vanishing pion mass of the 



masses of the fermionic and bosonic states in our present calculations are consistent with 
each other. Note that further away from the chiral limit the difference between fermionic 
and bosonic masses is expected to grow. 

In order to confirm that the remaining differences between the masses are finite lattice 
spacing effects, the calculations have to be extended to a third, even smaller lattice spacing, 
and extrapolated to the continuum limit. 

4 Conclusions 

In this work we present the current results of our simulations of TV = 1 supersymmetric 
Yang-Mills theory on a lattice. Our aim is to obtain a picture of the bound states of this 
theory in the nonperturbative regime. The results of our preparatory study on finite volume 
effects [S] are complemented with an analysis of the lattice artifacts. Both, finite volume 
effects and lattice artifacts increase the mass gap in the spectrum between the bosonic and 
the fermionic states. Hence the scales need to be chosen carefully for a reliable simulation 
of the theory. 

The parameters of our earlier results |3j have been chosen according to the experience of 
the QCD simulations. Our new, more detailed study has shown that these settings were on 
the safe side concerning the finite volume effects. The volume could even be reduced without 
a considerable systematic error. The supersymmetry breaking due to the discretisation 
effects or lattice artifacts, on the other hand, has been significant in our earlier results. The 
difference from the QCD expectations for the best simulation scales is not unexpected. The 
scale setting is ruled by the Compton wavelength of the particles under consideration. The 
volume has to be large and the lattice spacing has to be small in comparison to this length 
scale. In supersymmetric Yang-Mills theory there is no propagating particle corresponding 
to the lightest particle of QCD, the pion (the corresponding operator has been used only 
to tune the chiral symmetry restoration) . The absence of this light particle induces smaller 
finite volume effects, but obviously also larger discretisation errors. 

With our current parameters at the smaller lattice spacing the mass splitting is of the same 
order as the statistical and systematic uncertainties. For a quantitative analysis of lattice 
spacing effects and scaling an extrapolation to the continuum limit is required. In case of the 
gluino-glue this seems to be possible with the current accuracy of the results. For the other 
particles a reliable extrapolation requires further improvement of the measurements. The 
uncertainties of these measurements are in accordance with the expectations from QCD. 
The corresponding operators in QCD are also known to be hard to measure. In future 
studies we plan to improve the measurements and test the expected further reduction of the 
lattice artifacts. 
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